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I. INTRODUCTION 



Lattice QCD simulations with two light degenerate sea quarks (Nf = 2) as well as with a strange sea quark 
(Nf = 2 + 1) close to physical values of the pion mass are being carried out. Masses of low-lying hadrons are primary 
quantities that can be extracted using these simulations. Comparing the lattice and experimental values provides a 
check of lattice discretization effects. Such a comparison is necessary before one can use the lattice approach to study 
hadron structure. The European Twisted Mass Collaboration (ETMC) has generated a number of Nf = 2 ensembles 
at four values of the lattice spacing, ranging from 0.1 fm to about 0.05 fm, at several values of the light sea quark 
mass and for several physical volumes with maximally twisted mass fermions. We will use ensembles generated at the 
three smallest lattice spacings to evaluate the masses of strange and charm baryons. The strange and charm quarks 
are added as valence quarks. 

For heavy quarks the Compton wavelength of the associated heavy-light meson is comparable to presently attainable 
lattice spacings, which means that cut-off effects maybe large. The charm quark mass is at the upper limit of the range 
of masses that can be directly simulated at present. In order to obtain values for the masses that can be compared 
to experiment, it is important to assess the size of lattice artifacts. A first study of cut-off effects was carried out 
for light and strange baryons in Refs. 0, 0- I n this work we extend the study by including a finer lattice spacing 
and calculate besides the mass of strange baryons the masses of charm baryons. Having three lattice spacings the 
continuum extrapolation can be better assessed. 

In this work we compare our results in the strange baryon sector with recent results obtained with Clover-improved 
Wilson fermions with different levels of smearing. The PACS-CS Q and BMW Q collaborations evaluated the octet 
spectrum using two degenerate flavors of light quarks and a strange quark with mass tuned to its physical value. The 
PACS-CS has also computed the decuplet baryon masses. In addition, we compare with the LHPC that computed 
the octet and decuplet spectrum using a hybrid action with domain wall valence fermions on Kogut-Susskind sea 
quarks @ 

Besides the strange baryons, we also study the ground state spectrum of charm baryons with spin J = l/2 + and spin 
J = 3/2 + . Experimental searches of charm hadrons have received significant attention, mainly due to the experimental 
observation for candidates of the doubly charm baryons S+ (3520) and S+ + (3460) by the SELEX collaboration (^-Q. 
The 60 MeV mass difference between the singly and doubly charged states is difficult to understand since it is an 
order of magnitude larger compared to what is expected. No evidence was found for these states by the BABAR 
experiment [9[ and FOCUS Collaboration [Toj . The BELLE Collaboration [ll| finds S-states lower in mass, that 
can be candidates of excited states of S c but no doubly charm S. Additional experiments are planned at the new 
Beijing Spectrometer (BES-III) and at the antiProton ANnihilation at DArmstadt (PANDA) experiment at GSI, 
that can shed light on these charm baryon states. Several lattice QCD studies have been carried out to study charm 
baryons. We will compare the results of the current work with recent lattice QCD results all computed in a hybrid 
action approach where the charm valence quark was introduced on gauge configurations produced with staggered sea 
fermions by the MILC collaboration (l2HT^ . 




As in the case of the other lattice QCD studies of heavy baryons, also in this work we use a mixed action approach. 
For the strange and charm sector we use an Osterwalder-Seiler valence quark, following the approach employed in the 
study of the pseudo scalar meson decay constants [HI, [l6| . The bare strange and charm valence quark mass is tuned 
by requiring that the physical values of the mass of the kaon and D-meson are reproduced after the lattice results are 
extrapolated at the physical value of the pion mass. The ETMC Nf = 2 configurations [l?], [Hi analyzed in this work 
correspond to pion masses in the range of 260 to 450 MeV and three values of the lattice spacing corresponding to 
f3 = 3.9, 4.05 and 4.2 with r^/a = 5.22(2), 6.61(3) and 8.31(5), respectively. The Sommer parameter tq is determined 
from the force between two static quarks, the continuum value of which is determined to be 0.462(5) fm. At /3 = 4.2 
we use two ensembles, one corresponding to the lowest value of the pion mass considered in this work and one to 
the upper pion mass range. We find that the baryon masses, in general, show a very weak dependence on the lattice 
spacing and are fully compatible with an 0(a 2 ) behaviour with an almost vanishing coefficient of the a 2 term. This 
justifies neglecting the 0{a 2 ) term in extrapolating results to the continuum limit. 

An important issue raised by the twisted mass fcrmion formulation is isospin symmetry breaking. This symmetry, 
although exact in the continuum limit, is broken at a non-vanishing lattice spacing to 0(a 2 ). There are, however, 
theoretical arguments []j| and numerical evidences [2(J Hl[ that these isospin breaking effects are only sizable for the 
neutral pseudo scalar mass whereas for other quantities studied so far by ETMC they are compatible with zero. In 
this paper we demonstrate that also in the baryon sector these isospin breaking effects are in general small or even 
compatible with zero. Small isospin breaking effects decrease as the lattice spacing decreases and they vanish at the 
continuum limit. This corroborates our previous findings [H, HI]. The isospin breaking effects are relevant not only for 
neutral pions but also for other particles, e.g. the kaons. However, since the mass of the kaon is higher, the relative 
splitting (between K° and K + ) is less drastic. 

The paper is organized as follows: The details of our lattice formulation, namely those concerning the twisted mass 
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action, the parameters of the simulations, the interpolating fields used and the tuning of the strange and charm quark 
masses are given in Section II. Section III contains the numerical results of the baryon masses computed for different 
lattice volumes, lattice spacings and bare quark masses. Lattice artifacts, including finite volume and discretization 
errors, and continuum extrapolation are also discussed in Section III, with special emphasis on the 0(a 2 ) isospin 
breaking effects inherent to the twisted mass formulation of lattice QCD. The chiral extrapolations are analyzed in 
Section IV. Section V contains a comparison with other existing calculations. Our conclusions are finally drawn in 
Section VI. 



II. LATTICE FORMULATION 



A. The lattice action 



For the gauge fields we use the tree- level Symanzik improved gauge action [22| , which includes besides the plaquette 
term also rectangular (1 x 2) Wilson loops U\^ v 



^ = |E( 6 ° E {l-ReTr^Jl+fti Y, {l-ReTr(C7i^)} 



(1) 



jU,I/=l 



with b\ = —1/12 and the (proper) normalization condition bo = 1 — 8&i. Note that at b\ = this action becomes the 
usual Wilson plaquette gauge action. 

The fcrmionic action for two degenerate flavors of quarks in twisted mass QCD is given by 



S F = a A Yx{x){D w [U] + m + i^ 5 r 3 )x{x) 



(2) 



with t 3 the Pauli matrix acting in the isospin space, /i the bare twisted mass and the massless Wilson-Dirac operator 
given by 



1 



Dw[U] = -7 M (V 



where 



V M V(z) = - 



UUx)i/j(x + ap.) -ip(x) 



and V*V>(x) 



Ufj,(x — afi)ip(x — aji) — ip(x) 



(3) 



(4) 



Maximally twisted Wilson quarks are obtained by setting the untwisted quark mass too to its critical value m cr , while 
the twisted quark mass parameter \i is kept non-vanishing in order to work away from the chiral limit. In Eq. ^ the 
quark fields x are m the so-called " twisted basis" . The " physical basis" is obtained for maximal twist by the simple 
transformation 



(5) 



ip(x) = exp [^75T J x(x), ip(x) = x(x) exp ( — j 5 t 
In terms of the physical fields the action is given by 

sf = « 4 £^) (^7m[v m + vy - *7 5 r 3 (-y v M v; 



H ip{x) 



(6) 



In this paper, unless otherwise stated, the quark fields will be understood as "physical fields", ip, in particular when 
we define the baryonic interpolating fields. 

A crucial advantage of the twisted mass formulation is the fact that, by tuning the bare untwisted quark mass 
too to its critical value m cr , all physical obscrvables arc automatically 0{a) improved. In practice, we implement 
maximal twist of Wilson quarks by tuning to zero the bare untwisted current quark mass, commonly called PCAC 
mass, topcac [H|j which is proportional to mo — m cr up to 0(a) corrections. The value of m cr is determined at each 
P value at the lowest twisted mass used in our simulations, a procedure that preserves 0(a) improvement and keeps 
0(a 2 ) small [l8l |2~H]. The twisted mass fermionic action breaks parity and isospin at non- vanishing lattice spacing, as 
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it is apparent from the form of the Wilson term in Eq. ([5]) . In particular, the isospin breaking in physical observables 
is a cut-off effect of 0(a 2 ) (25[. To simulate the strange quark in the valence sector several choices are possible. 
The strange and charm quarks are added as Osterwaldcr-Sciler valence quarks and their action reads 

c 

S£? avy = « 4 E E ^ ( y( V '< + V P - i V ^M + M or + *75 Mft) Xh(x) (7) 
x h=s 

where [i s and /i c are the strange and charm valence quark masses. This is naturally realized in the twisted mass 
approach by introducing two additional doublets of strange and charm quarks and keeping only the positive diagonal 
component of T3. The mo value is taken to be equal to the critical mass determined in the light sector, thus guaranteeing 
the O(a) improvement in any observable. The reader interested in the advantage of this mixed action in the mesonic 
sector is referred to the Refs QJ, [U, f^MM- 



B. Simulation details 



The input parameters of the calculation, namely /?, L/a and an are summarized in Table HI The corresponding 
lattice spacing a and the pion mass values, spanning a mass range from 260 to 450 MeV, are taken from Ref. j29j . 
At w 300 MeV we have simulations for lattices of spatial size L = 2.1 frn and L = 2.7 fm at f3 = 3.9 allowing 
to investigate finite size effects. Finite lattice spacing effects are investigated using three sets of results at f3 = 3.9, 
f3 = 4.05 and (3 = 4.2. These sets of gauge ensembles allow us to estimate all the systematic errors in order to have 
reliable predictions for the baryon spectrum. 







P = 4.2, 


a = 0.056(1) fm r /a = 8.31(5) 




32 3 


x 64, L 


= 1.8 fm 


a^ sca 0.0065 
Statistics 240, 76 

(GeV) 0.4698(18) 
m n L 4.24 




48 3 


x 92, L 


= 2.7 fm 


a/Xsea 0.0020 
Statistics 458, 456 
(GeV) 0.262(1) 
mn^L 3.55 








P = 4.05, 


a = 0.070(1) fm, r /a = 6.61(3) 




32 3 


x 64, L 


= 2.13 fm 


ap a ea 0.0030 0.0060 0.0080 
Statistics 144, 144 194, 193 201, 201 
m n (GeV) 0.2925(18) 0.4035(18) 0.4653(15) 
m n L 3.31 4.57 5.27 








P = 


= 3.9, a = 0.089(1) fm, r /a = 5.22(2) 




24 3 


x 48, L 


= 2.05 fm 


a/i S ea 0.0040 0.0064 0.0085 
Statistics 4112, 310 545, 278 1817, 369 


0.010 

477, 475 








(GeV) 0.3032(16) 0.3770(9) 0.4319(12) 0.4675(12) 
m^L 3.25 4.05 4.63 5.03 


32 3 


x 64, L 


= 2.74 fm 


a^ SC a 0.0030 0.0040 
Statistics 659, NA 242, NA 
(GeV) 0.2600(9) 0.2978(6) 
m^L 3.74 4.28 





TABLE I: Input parameters (/3,L,/i) of our lattice simulations and corresponding lattice spacing (a) and pion mass (m^). The 
statistics refer to the number of configurations used in calculation of the masses of the strange and charm baryons. The first 
entry gives the number used for the tuned value of the strange quark and the second for the tuned value of the charm. An 
entry NA indicates that no masses were computed. The lattice spacing was determined using the nucleon mass [2!|. 



C. Tuning of the bare strange and charm quark masses 

The de pen dence of the pseudoscalar meson mass on the valence and sea quarks can be written as a polynomial of 
the form |l5j 

a 2 Mp S (a^ sca , a/ii, ^ 2 ) = B (afii + a(j, 2 )[l + a v £,i2 + a sca ^ sca + avvth + a sca£sca + av,sca£,i2^sca + ayuCi)^] > ( 8 ) 
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where /i soa is the sea quark mass and am^ a m are the valence quark masses, £3 = B^aml '(4tt f) 2 , £y = 2B^{am + 
a/j,j)/(4:Trf) 2 and £oy = Bo(am — afj,j)/(4:irf) 2 . For the /? = 3.9 ensembles we consider in total 164 pseudoscalar 
meson masses using all possible combinations of sea and valence quark masses. Namely, we have considered 150 
combinations obtained from a^t sea and am independently taking the values 

{a^ sca ,a^i} = {0.0040 0.0064 0.0085 0.0100 0.0150} 

whereas am takes the values 

am = {0.0040 0.0064 0.0085 0.0100 0.0150 0.0220 0.0270 0.0320}. 
We have an additional 12 combinations coming from the combinations 

a^sca = Mi = {0.0040 0.0064 0.0085 0.0100} , m = {0.24 0.27 0.30}, 
plus two extra combinations from 

a^ sca = ani = 0.0040 , am = {0.0217 0.25} . 

For the tuning at j3 = 4.05 we use the following 20 combinations 

a^sca = am = 0.0030 am = {0.0030 0.014 0.0166 0.020 0.17 0.20 0.23 0.26} 

a^sca = am = 0.0060 am = {0.0060 0.0166 0.019 0.025} 

a^ sca = am = 0.0080 am = {0.0080 0.014 0.0166 0.020 0.17 0.20 0.23 0.26}. 
For the tuning at j3 = 4.2 we consider 10 pseudoscalar meson masses: 

a// sca = am = 0.0065 am = {0.0065 0.14 0.16 0.185 0.21} 

aMsca = am = 0.0020 am = {0.0020 0.012 0.015 0.136 0.17}. 

In Figs. [1] and [21 we show representative fits to the pseudo-scalar masses in the range of the kaon and D-meson 
masses using the expression given in Eq. (|8"|). The values of the strange and charm quark masses are varied until the 
resulting kaon and D-mcson masses are matched to their physical values. The resulting fit parameters are listed in 
Table [Til We note that for (3 = 3.9 two fitting ranges are used, one range spanning the the strange quark mass and 
one the charm quark mass, For j3 = 4.05 and (3 = 4.2 we fit all data together since we do not have enough mass 
combinations in order to apply Eq. (J8J). If one does the same for j3 = 3.9 then the tuned value for the strange quark 
mass is a/j, s = 0.0216(7) compatible with the value of afi s = 0.0217(5) if we restrict the fit to the strange region. In 
addition, at each /3-value we can restrict the fit in the charm region using the Ansatz 

m B = m% + bm + c/m- (9) 

The tuned charm quark value is found to be compatible with the one extracted using Eq. ([5]) . This procedure can be 
carried out using either the lattice spacing determined from the nucleon mass or form f„ . The difference in the tuned 
masses reflects the systematic error in setting the scale. 

The tuned values of the strange and charm quark masses m s and m c obtained at the the physical pion are given in 
Table HTT1 In a previous paper, the ETMC computed pseudo-scalar meson masses for a number of sea and valence quark 
masses using the j3 = 3.9 gauge configurations. Matching the experimental value of the mass ratio of the kaon to the 
pion, rrifc jra^ , the bare strange quark mass was determined (l5| . Depending on the polynomial fit used the values for 
am s at j3 = 3.9 varied from 0.0243(5) to 0.0218(10). Thus, our value of a/j, s = 0.0216(7) from matching the physical 
value of the kaon mass in combination with the lattice spacing determined from the nucleon mass is compatible with 
the value determined in Ref. [15j. Such an agreement is satisfactory and shows that the two procedures lead to the 
same determination within the uncertainties associated with the extrapolation. The systematic error introduced from 
the way the lattice scale is fixed can be assessed by comparing the tuned values extracted using the lattice spacing 
determined from the nucleon mass and from the pion decay constant f n . The values of the lattice spacing determined 
using f n , taken from Ref. [H, are 0,3=3.9 = 0.0801(14) fm, a p=i . 05 = 0.0638(10) fm and a p=4 . 2 = 0.05142(38) fm. In 
Table [TTT1 we give the tuned values for the charm and strange quark masses expressed in physical units. As can be 
seen, the values for the charm quark masses are in agreement, whereas for strange quark masses the differences are 
within about two standard deviations. 
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p = 3.9 (strange quark) 


p = 3.9 (charm quark) 


P — 4.05 


P = 4.2 




2.zo2(o) 


1. 38(D) 


l.D02(0 ) 


1.290(0 ) 


f 


0.077(2) 


0.112(2) 


0.093(4) 


0.069(4) 


av 


-0.45(2) 


0.3(1) 


0.85(5) 


0.56(5) 


Q-sea, 


0.0 


0.0 


0.0 


0.0 


avv 


3.0(1) 


1.8(5) 


-4.0(3) 


2.6(2) 


& V.sea, 


0.0 


0.0 


0.0 


0.0 


^sca 


0.0 


0.0 


0.0 


0.0 


ClVD 


-2.25(3) 


-1.4(6) 


4.94(5) 


-1.8(3) 


xVd-of 


0.51 


1.33 


4.52 


4.40 



TABLE II: The values of the fit parameters. 
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=0.003 
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=0.008 
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FIG. 1: Pseudo-scalar meson masses for P = 4.05 as a function of the heavy valence quark mass afi2 in the relevant mass range 
for the strange quark (left) and charm quark (right). In all the examples shown the sea quark mass a/i SC a is set equal to the 
light valence quark mass a/ii. 




FIG. 2: The dependence of the pseudo-scalar masses on the strange quark mass at P = 3.9 (left) and charm quark mass at 
P — 4.05 (right). The solid line shows the variation of the pseudoscalar mass at the tuned strange (left) and charm (right) 
quark masses. The dashed vertical line corresponds to the value of afii at which the physical pion mass is recovered. The 
asterisk denotes the physical value of the kaon and D-meson mass in lattice units. The lattice spacing is determined from the 
nucleon mass. 



D. Interpolating fields 

The low-lying baryons belonging to the octet and decuplet representations of SU(3) are given in Figs. [3] and [4] 
respectively. They are classified by giving the isospin, /, the third component of the isospin, 7 3 , the strangeness (S), 
the spin and the parity. In order to extract their masses in lattice QCD, we evaluate two-point correlators. We use 
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p 






H* (GeV) 


Mc (GeV) 


l*i*(GeV) 


I& (GeV) 


3.9 


0.0216(7) 


0.27(3) 


0.0478(16) 


0.598(66) 


0.0431(17) 


0.64(12) 


4.05 


0.0178(5) 


0.21(1) 


0.0501(14) 


0.591(28) 


0.0451(12) 


0.556(31) 


4.2 


0.014(1) 


0.17(2) 


0.0493(35) 


0.598(70) 


0.0464(15) 


0.575(38) 



TABLE III: The strange and charm quark masses at each value of /3 tuned using the kaon and D-meson masses and the lattice 
spacing determined from the nucleon mass are given in lattice units in the second and third columns respectively. The tuned 
strange and charm quark masses in GeV are given in the fourth and fifth columns. In the sixth and seventh columns we 
determined the corresponding masses in GeV using the lattice determined from /„■. 



interpolating fields to create these states from the vacuum that have the correct quantum numbers and reduce to the 
quark model wave functions in the non-relativistic limit. The interpolating fields used in this work are collected in 
Tables 23 and [V] [HI, for the octet and decuplet respectively. 

Charm baryons with no strange quarks are obtained from the interpolating fields of strange baryons by replacing 
the strange with the charm quark. There are additional charm baryons containing strange quarks, giving a 20-plet 
of spin-1/2 and a 20-plct of spin-3/2. In most of this work we do not consider the particles that contained both a 
strange and charm quarks. For the lattice with the smallest lattice spacing and at the smallest pion mass we also 
consider the spin-1/2 S c , E' c , r2 c and fl cc and the spin-3/2 S*, 0* and £2* c . The interpolating fields for these baryons 
are given in Table IVII 



s = 



= -1 E" 



s = -2 




q = l 



s = Q 



s = -1 



s = -2 



s = -3 



A+ A+^ 




q = 2 



9 = 1 



q = Q 



q = -l 



FIG. 3: The low lying baryons belonging to the octet rep- 
resentation labeled by the value of 73 and hyper-charge. 



FIG. 4: The low lying baryons belonging to the decu- 
plet representation labeled by the value of Is and hyper- 
charge. 



Strangeness 


Baryon 


Interpolating field 


I h 


5 = 


P 
n 


X U = ^abc(d,lC^5Ub)d c 


1/2 +1/2 
1/2 -1/2 




A 


X AS = -^^abc{2(ul C^sdb)s c + {u^C^sSb)d c - (dlC-y 5 Sb)u c } 





5 = 1 


£+ 
E° 


X E = £abc{ulC'y'oSb)Uc 

X E ° = -^^abc{(ulC^ 5 s b )d c + (dlC"/ 5 s b )u c } 


1 +1 

1 +0 




sr 


X E = £abc(da Cj5Sb)d c 


1 -1 


5 = 2 


w0 


X" — £abc(S a CysUbjSc 
X" = £abc(SaC-y5db)s c 


1/2 +1/2 
1/2 -1/2 



TABLE IV: Interpolating fields and quantum numbers for the baryons in the octet representation. 



Local interpolating fields are not optimal for suppressing excited state contributions. We instead apply Gaussian 
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Strangeness 


Baryon 


Interpolating field 


I 


I z 




A++ 


A++ / Try \ 
Xp, = ^abciUaCj^UbjUc 


3/2 


+3/2 


S = 


A+ 


Xm + = -^£abc{2(iiaC'ypd b )uc + (m„ Cj fl Ub)d c } 


3/2 


+ 1/2 




A 

A u 


Xm = 75 £ a&c{2(d a Cf^Ubjdc + (d a C"fp.d b )u c ) 


3/2 


-1/2 




A" 


Xm = e a bc(dlC^d b )d c 


3/2 


-3/2 






Xm = sJle abc {{u Ta C ltl u h )s c + 2(u Ta C 1 pS b )u c } 


1 


+ 1 


S=l 


E*° 


X f ° = y^e o6c {(u To C7 # ,dV + (d^Oy^V + (s^C^V } 


1 


+0 




E*~ 


Xf" = A /fe a,,c {(d Ta C7 M d i, ) S c + 2(d T[l C7 MS i ')d c } 


1 


-1 


S = 2 




Xm = eabc(s a C7 M Ui,)s c 


1/2 


+ 1/2 






Xm = ea& c (srC7 M di,)sc 


1/2 


-1/2 


S = 3 


Q- 


Xm = tabdsIC-y^s^Sc 





+0 



TABLE V: Interpolating fields and quantum numbers for baryons in the decuplet representation. 



J = 1/2 


J = 3/2 


X" c = -^eabc{2(SaC'ysd b )cc + (slCjs^dc - (d£ C-ysc b )s c } 


Xm c = 


^e abc {{s Ta C^d b )c c + (d Ta C lllC »)s c + (c Ta C^s»)d c } 


X~" = -^^aba{(SaCj 5 Cb)d c + {dlC~l 5 C b )s c } 






X^ = ^abc{slC-y 5 C b )s c 




X f = Jle^{(s T °Cj,s b y + 2(s T «C>y,c b )sC} 


X " = ^abc{clC^ 5 Sb)c c 




Xm" = CabciCaClpS^Cc 



TABLE VI: Interpolating fields for the spin-1/2 S c , Hj,, f2 c and fi cc , and the spin-3/2 H*, Q* and $l* c baryons. 



smearing to each quark field, q(x, t): g smear (x, t) = J2 y ^( x i Yi U(t))q(y, t) using the gauge invariant smearing function 

F(x, y; C/(t)) = (1 + off) n (x, y; £/(*)), (10) 

constructed from the hopping matrix, 

H (x, y; tf(t)) = ^ (V,(x, tj^y-i + E/?(x - i, t)S x . y+l ) . (11) 
i=i ^ ' 

Furthermore we apply APE smearing to the spatial links that enter the hopping matrix. The parameters of the 
Gaussian and APE smearing are the same as those used in our previous work devoted to the nucleon and A masses [l| . 



E. Two-point correlators 

To extract masses in the rest frame we consider two-point correlators defined by 

C±{t,p = 0) = -Tr(l ± 7o ) {Jx(xdnk,tank)J X ( ^source ) ^source)) j ^ ^sink ^source • (1^) 

x sink 

Space-time reflection symmetries of the action and the anti-periodic boundary conditions in the temporal direction 
for the quark fields imply, for zero three-momentum correlators, that C~t (t) = —C^(T — t). So, In order to decrease 
errors we average correlators in the forward and backward direction and define: 

C x (t) = C+(t)-C x (T-t). (13) 

In order to decrease correlation between measurement, we choose the source location randomly on the whole lattice 
for each configuration. Masses are then extracted from the so called effective mass which is defined by 

mfgft) = -log(C x (t)/C x (t-l)) = m x +log ( ^^(t-t) J > ( 14 ) 
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where A j = ?ti,; — mx is the mass difference of the excited state i with respect to the ground mass mx ■ 
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FIG. 5: Representative effective mass plots for j3 = 4.2 and afii = 0.002. For the strange baryons (Left) we used afi a = 0.015, 
while for charm baryons (Right) a/i c = 0.17. 

In Fig. [5]we show representative examples of the effective masses of strange and charm baryons. As can be seen, a 
plateau region can be clearly identified. What is shown in these figures are effective masses extracted from correlators 
where smearing is applied both at the sink and source. Although local correlators are expected to have the same value 
in the large time limit, smearing suppresses excited state contributions yielding a plateau at earlier time separations 
and to a better accuracy in the mass extraction. We therefore extract the masses using smeared source and sink. Our 
fitting procedure to extract mx is as follows: The sum over excited states in the effective mass given in Eq. ([14"]) is 
truncated keeping only the first excited state. Allowing a couple of time slice separation the effective mass is fitted 
to the form given in Eq. (|14j) . This yields an estimate for the parameters c\ and Ai. The lower fit range is increased 
until the contribution due to the first excited state is less than 50% of the statistical error of mx ■ This criterion is 
in most of the cases in agreement with a x 2 /d.o.f. < 1. In the cases in which this criterion is not satisfied a careful 
examination of the effective mass is made to ensure that the fit range is in the plateau region. 



III. LATTICE RESULTS 



Before we extrapolate our lattice results on the strange and charm baryon masses to the physical point, we need to 
examine their dependence on the heavy quark mass as well as cut-off effects. We collect lattice results for the masses 



of the strange and charm baryons in the Appendix. The errors are evaluated using jackknifc and the T-method 34 1 
to check consistency. 

In Figs. [5] and [7] we show the dependence of the strange and charm baryon masses on the the strange and charm 
quark mass, respectively. Overall, the data display a linear dependence on both the strange and charm quark mass. 
One can therefore interpolate between different values of quark masses, if needed. 
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FIG. 6: Baryon masses versus the strange quark mass for f3 = 3.9 and a/i; = 0.0064 (circles), j3 = 4.05 and afii — 0.006 
(squares), /3 = 4.2 and a/i; = 0.0065 (triangles). 



A. Strange baryon mass with strange quark mass tuned to its physical value 

In this section we restrict our analysis only to the subset of data obtained at the tuned values of the strange quark 
mass. Namely, for j3 = 3.9 and j3 = 4.05, we use the tuned value given in Table IIIIl whereas for j3 = 4.2 we use 
afi s = 0.015 which agrees with the tuned strange quark mass within error bars. 

It is interesting to examine the degree of isospin splitting as a function of the lattice spacing. The splitting is 
expected to be zero in the continuum limit. In Fig. [5] we show the mass of S + , S° and E~ at ft = 3.9,4.05 and 4.2. 
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As expected the mass splitting among the three charge states of the £ decreases with the lattice spacing. The same 
behavior is observed for the other strange particles studied in this work. This is shown in Fig. [S] where we plot the 
mass difference as a function of a 2 at our smallest and heaviest pion mass. As can be seen, the mass difference is 
consistent with zero for all particles at the smallest lattice spacing. The small non-zero values seen for the £ and 5 
particles are just outside one standard deviation. Therefore, the general conclusion is that indeed isospin splitting is 
small at these values of the lattice spacing and it vanishes at the continuum limit. Since for finite a there are small 
differnccs, for the chiral extrapolation where we use all lattice data we do not average the masses for the different 
charge states of the £, S and S*. 
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Volume effects can be studied at j3 — 3.9 where we have simulations at two volumes at pion mass of about 300 McV. 
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FIG. 9: Isospin mass splittings in the strange-quark sector at the smallest (left) and largest (right) pion masses used in our 
computation. 
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As can be seen in Figs. 151 ITUl and [TT1 results at different volumes are consistent. Therefore, we conclude any volume 
effects are smaller than our statistical accuracy. 
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FIG. 10: Left: The strange octet baryon mass in units of ro versus (a/ro) 2 for several particles in the octet and decuplet at 
two hxed values of m^ro. The solid line is a linear fit in (a/ro) 2 and the dashed line to a constant. Right: Chiral extrapolation 
at fixed strange quark mass. Results at /3 = 3.9 are shown with the open circles for L = 2.1 fm and with the filled circles for 
L = 2.7 fm, at j3 = 4.05 with the open squares and at /3 = 4.2 with the open triangles. 



In order to examine the continuum limit we interpolate our lattice result at a given pion mass in units of ro. 
For /3 = 4.2 we have simulations at only two values of the pion mass at the upper and lower range of pion masses 
considered in this work, namely for = 0.254 GeV and m w = 0.459 GeV. Therefore we interpolate the results for 
the other two values of j3 to these two pion masses. In Figs. [TU] and [TT1 we show results for the octet and decuplet 
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FIG. 11: Left: The strange decuplet baryon mass in units of ro versus (a/ro) 2 at two fixed values of m^ro. The solid line is 
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of the E* we show the three charge states E* + (purple), E*° (blue), and E*~ (red). The spread of data at a given pion mass 
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strange baryon masses, respectively, at our three values of the lattice spacing. We perform a continuum extrapolation 
by perfroming a linear fit in (a/ro) 2 as well as to constant. As can be seen from Fig. [101 the values obtained in the 
continuum limit agree for all octet baryons. In the case of the decuplet the statistical errors are larger and the value 
obtained at a = with the linear fit carries a large error. The value obtained using a constant fit has a smaller error 
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and it is compatible with one the one obtained using a linear fit. Therefore, for a given charge state and within the 
current statistical accuracy, the (a/ro) term can be taken as negligible. Therefore, we can use results at all /3-values 
to extrapolate to the physical point since cut-off are small for a given charge state. There are two exceptions in the 
case of the decuplet. At ft = 4.05 the mass of the S* and O at the lowest pion mass are systematically higher than 
that at the other two ft- values. Since the results at ft = 3.9, with larger lattice space are consistent with those at 
ft = 4.2 we conclude that this is not a cut-off effect. 



B. Charm baryon mass with charm quark mass tuned to its physical value 

As in the previous subsection we consider results obtained at the tuned charm mass given in Table IIIII The 
only exception is at ft = 4.2 at the heavy pion mass, where we have results close to the tuned value, namely at 
a/i c = 0.16, 0.185, 0.21. As we have seen, the dependence on the heavy quark is linear and therefore the charm baryon 
masses at the tuned value can be easily determined by a linear interpolation. 

We follow the same analysis as in the case of the strange baryon sector. In Fig. rjj] we show the mass difference 
between different charged states as a function of the lattice spacing at the smallest and largest pion masses used in 
this work. As can be seen, the mass splittings are zero at the smallest value of the lattice spacing for all particles 
confirming restoration of isospin symmetry in the continuum limit. Furthermore, except for the case of the S cc mass, 
the mass splitting is consistent with zero also at the other two /3-values. Therefore, for all particles except the S cc 
one may average over the mass of different charge states. 

In order to examine the continuum limit we interpolate our results at the three /3-values at a given pion mass in 
units of ro- In Figs. fiUl and [T4l we show the mass in the octet and decuplet charm sector as a function of lattice 
spacing for a given charge state, at the smallest and largest value of the pion mass. A linear fit in (a/ro) 2 and a 
constant fit yield consistent results at the continuum limit, albeit with large errors in the case of th elinear fit. We 
also note that at the largest pion mass, although results at ft = 3.9 are in agreement with those at ft = 4.2 indicating 
neglegible C(a 2 )-dependence, at ft = 4.05 the results are systematically below. We note that we show only statistical 
errors. Systematic errors due, for example, to the matching are not shown. As discussed in the next section these 
are (5-10)%. Therefore, a reasonable way to extrapolate our results in the charm sector is to compare the chiral 
extrapolation using all lattice data to those using results at ft = 3.9 and ft = 4.2. We will take the different between 
the two values at the physical point as an estimate of a systematic error. 

IV. CHIRAL EXTRAPOLATION 

Having determined that 0(a 2 ) effects are small for the lattice spacings considered here we can combine our lattice 
results at the various /3-values to extrapolate to the physical pion mass (physical point). 

For the strange baryon sector, we consider SU(2) heavy baryon chiral perturbation, which was found to describe 
lattice data satisfactorily Q . To leading one- loop one can described the pion mass dependence using 

tub = tUb — 4c^TO 2 r + cm% (15) 

where c is a known coefficient given in Ref. [2[ . For completeness we give below the coefficients c [H, [36| . For the 
octet baryons A, £ and S: 

ffL 2. 9 | s+g | s /3 3ffj s 

167T/2 ' 167T/2 ' 167T/2 

respectively, and for the decuplet baryons E*, S* and f2: 

_ 10 5 |» E » 5 5 l. a . 



9 16./- 316,/^-°' 

In addition we consider next to leading order SU(2) xPT results (37j- The expressions are included here for 
completeness: 



2 a 2 

4 L °{m.) = mf-^ml _ ^( m „ A AS , A) - ^(m., A AS . , A) 



16tt/2 (4tt/ t ) 2 
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FIG. 12: Isospin mass splitting in the charm sector at the smallest and largest values of the pion mass used in this work. 
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FIG. 13: Left: The mass of octet charm baryons in units of ro versus (a/ro) 2 for several particles in the octet and decuplet at 
two fixed values of the pion mass m n ro. Right: Chiral fits to data at f3 — 3.9 (circles) and at /3 — 4.2 (triangle) are shown with 
the green band and to all data with the lue band. For the H cc we have extrapolated whereas for E c we averaged over E+ + 
(purple), E+ (blue) and E° (red). The spread of data at a given pion mass with the same symbol indicates the level of isospin 
breaking for E. 
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FIG. 14: Left: The mass of decuplet charm baryons in units of ro versus (a/ro) 2 at two fixed values of m n ro. Right: Chiral fits 
to the data at ft = 3.9 (circles) and at j3 = 4.2 (triangle) are shown with the green band and to all data with the blue band. For 
H* c we averaged over 3*^ + (blue) and H*,t (red) and for E* we averaged over (purple), E* + (blue) and E*° (red). The 

spread of data at a given pion mass with the same symbol indicates the level of isospin breaking for the and H* c particles. 
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Particle(PDG) 


(GeV) 


-4c^ (GeV" 1 ) 


X 2 /d.o.f. 


m (GeV) 


£"(1193) 


1.1368(70) 


3.560(40) 


2.7 


1.1930(62)(660) 


=^(1315) 


1.3334(46) 


1.386(26) 


0.82 


1.3538(41)(179) 


A(1116) 


1.0678(64) 


4.362(37) 


1.04 


1.1276(57) (721) 


EL (1384) 


1.4244(58) 


2.807(34) 


2.4 


1.4757(51) (740) 


H*"(1531) 


1.4808(96) 


1.582(58) 


3.3 


1.5113(89)(400) 


f2(1673) 


1.7522(76) 


0.361(45) 


2.0 


1.7591(67) (200) 



TABLE VII: The bare mass and cb (related to the a-term by as = — 4csm^.) determined from fitting to the NLO expressions 
for strange baryons at the tuned strange quark mass. In the last column, we give the mass in GeV that we obtain at the 
physical point using the NLO expressions. The second error given in the parenthesis is an estimate of the systematic error 
coming from a comparison between the values obtained at the physical point using the leading order (LO) expressions given in 
Eqs. 1 161 and 1 171 and the NLO given by Eqs. 1181 and ll9l In the case of the Q the systematic error is estimated by evaluating the 
impact of the error of the tuned strange quark mass on the extrapolated Q mass. 



and for the decuplct baryons: 



mi, L °{m«) = m^-iSml-l^m^-j^^ ^(m„, -A EE ., A ) 

mg LO (m w ) = m^~4c^ml (19) 
with the non analytic function f38j ] 



» ,x , o /To n , /A - VA 2 - to 2 + ie \ 3 A fm 2 \ - /4A 2 \ . . 

' A ) = ^ " AVA 2 -™ 2 + ie log [ A + VA2 ^ m2 + ie ) - Am 2 log J - A^ log J (20) 

depending on the threshold parameter Axy = — m j? an( f 011 ^ ne sca l e ^ of chiral perturbation theory, fixed to 
A = 1 GeV. For A > the real part of the function F(m, A, A) has the property 



-T{m, A, A) m < A 

-F(m, A, A) + 2tt (m 2 - A 2 ) 3/2 m > A 



7-(m,-A,A) = <j ^LV^.n/a a2\^/ 2 ^ a ( 21 ) 



which corrects a typo in the sign of the second term in Ref. @. We follow the procedure of Ref. Q and fix the 
nuclcon axial charge qa and pion decay constant / ff to their experimental values (we use the convention such that 
frr = 130.70 MeV). The remaining pion-baryon axial coupling constants are taken from SU(3) relations [37| ■ The fit 
parameters extracted for fitting to the NLO are given in Table I VIII The deviation of the mean values obtained at the 
physical point when the results are fitted to leading order i.e. to Eq.( [T5|) with c = and when they are fitted to the 
NLO expressions provide an estimate of the systematic error due to the chiral extrapolation. We give this error in 
Table IVHI In the case of the il there is no difference between leading order and next to leading order. Since the f2 
contains three strange quarks any systematic error in the tuning of the strange quark mass will be largest in this case. 
Having results at several values of the strange quark mass we can estimate the change in the fl mass if the strange 
quark mass takes the maximum and minimum value allowed by the statistical error in the tuned strange quark mass. 
We take the difference in the mean values at the physical point obtained by varying the strange quark mas to be the 
systematic error due to the tuning. In Table IVTI1 we give the systematic error on the mass of SI that we find following 
this procedure. This gives an upper bound of the error expected from the uncertainty in the tuning. As can be seen 
this is smaller as compared to the systematic error due to the chiral extrapolation and therefore it is only taken into 
account for the case of the fl. 

In Figs. [TU1 and ITT1 we show the chiral extrapolation for the octet and the decuplet. In the case of the S and S* the 
physical point is reproduced. However, for most other particles the lattice results extrapolate to a higher value. The 
worse deviation is seen for the f2. Since this has three strange quarks it may indicate that the tuning of the strange 
quark mass performed using the kaon mass introduces a systematic error. One can study partial quenching effects 
using twisted mass fermion simulations with a dynamical strange quark. This will be considered in a future study. 

In the charm baryon sector we use the Ansatz 

m B = a + bml + cm*, (22) 
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Particle(PDG) 


m (GeV) 


Am (GcV) 


E C , Q „ (2.454) 


2.494(47) 


0.143 


1 — 'CC 


3.563(25) 


0.397 


AJ(2286) 


2.229(43) 


0.223 


Xtav (2.520) 


2.650(39) 


0.147 


1 — 'CC,CMJ 


3.672(42) 


0.274 




4.702(11) 


0.308 



TABLE VIII: For each particle listed in the first column we give in the second column its mass at the physical pion mass using 
for the chiral extrpoaltion the masses computed at the tuned value of the charm quark mass m c . In the third column we give 
the mass difference between the baryon masses obtained at the tuned value of m c and at the tuned value plus the error, after 
extrapolation to the physical point. This is done at f3 = 3.9 where we have computed the masses at m c ± error. 



Particle(PDG) 


m u B (GeV) 


-4c s (GeV" 1 ) 


c (GeV"") 


rVd.o.f. 


m (GeV) 


E c , a „ (2.454) 


2.437(25) 


1.92(54) 


-2.09(91) 


1.1 


2.468(17)(23) 


1 — 'CC 


3.476(35) 


2.39(83) 


-3.39(1.5) 


2.7 


3.513(23)(14) 


AJ(2286) 


2.198(40) 


2.99(96) 


-3.6(1.7) 


0.10 


2.246(27)(15) 


E*, a „ (2.520) 


2.520(25) 


2.37(51) 


-2.96(86) 


1.3 


2.556(18)(51) 


1 — 'cc,av 


3.571(25) 


2.02(57) 


-2.62(99) 


1.0 


3.603(17)(21) 




4.6706(53) 


0.327(35) 


0. 


2.5 


4.6769(46) (30) 



TABLE IX: Parameters of the chiral fit for charm baryons at the tuned charm quark mass fitting results at /3 = 3.9 and /3 = 4.2. 
The last column is our prediction (in GeV) at the physical point. The statistical error is given in the first paranthesis and the 
systematic, computed by comparing the fit with all lattice data, in the second parenthesis. 

motivated by SU(2) HBxPT to leading one-loop order, with c taken as a fit parameter. For the f2 ccc we set c = 
since one does not expect a cubic term. 

In order to assess the systematic error associated with the tuning of the charm quark mass wc consider our results 
at [3 = 3.9. At this value of (3 we have computed the charm baryon masses at the tuned value of the charm quark 
and at values of the charm quark shifted by the error on the tuned value. Since these computations were performed 
at four different light quark masses we can perform a chiral extrapolion using using the Ansatz of Eq. (|22l for the set 
of masses obtained at the tuned value and at the value shifted by the error. The difference in the masses obtained at 
the physical pion mass is given in Tab. IVIIII As can be seen, this difference intorduces an error that varies between 
about 5% and 10%. This gives an estimate of the systemastic error due to the tuning of the charm quark mass. Since 
this analysis can only be done at /3 = 3.9 we can only make a qualitative estimate of this error. Therefore in what 
follows we will not quote this error on our values,. However, one has to bear in mind that our final values can have a 
systematic error of about 10% due to the tunning. 

In Figs. [13] and [14] we show fits at for our three /3 values. We show fits using all data and fits using only data at 
f3 = 3.9 and j3 = 4.2. The latter case yields a better fit with a smaller value of x/d.o.f. and this is the value quoted in 
Table IIXI This is particularly noticeable for the case of il ccc where the results at (3 = 4.05 are systematically lower. 
This maybe due to a small mismatch in the tuned value of the charm quark mass, which for the fl C cc that contains 
three charm quarks would lead to the largest deviation. We take the difference in the extrapolated values at the 
physical points when we exclude the f3 = 4.05 data from the fit as a systematic error. 

The extrapolation of the lattice data reproduce the mass of experimentally measured charm baryon masses within 
a standard deviation, namely the mass of the E c , the A c and the £*. Therefore, the extrapolated lattice value can be 
taken as a prediction for the mass of the S* c and the fi C cc, within one standard deviation. 

V. COMPARISON WITH THE RESULTS OF OTHER LATTICE FORMULATIONS 

In this section we compare our results with those using different discretization schemes by other collaborations. We 
also include a comparison for the nucleon and A masses although they were not discussed in detail until now. 

Several collaborations have calculated the strange baryon spectrum. The Budapest-Marseille- Wuppertal (BMW) 
collaboration carried out simulations using tree level improved 6-step stout smeared Nf = 2 + 1 clover fermions and 
a tree level Symanzik improved gauge action. Volume effects were studied using lattices of spatial extent of 2 fm to 
4.1 fm. The continuum limit was taken using results produced on three lattice spacings of a = 0.065 fm, 0.085 fm and 
a = 0.125 fm. Using pion masses down to 190 MeV a polynomial was performed to extrapolate to the physical value 
of the pion mass The PACS-CS collaboration obtained results using Nf = 2 + 1 non-perturbatively 0(a) improved 
clover fermions on an Iwasaki gauge action on a lattice of spatial length of 2.9 fm and lattice spacing a = 0.09 fm 
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FIG. 15: The results of this work for the octet strange baryons are shown with the filled (red) circles, results using clover 
fermions by the BMW Collaboration are shown with the open (blue) square 0] and by the PACS-CS Collaboration with 
the filled (magenta) triangles [3], and domain wall valence quarks on a staggered sea by the LHPC with the open (green) 
triangles [5(. The experimental value is shown with the asterisk. 



The QCDSF-UKQCD collaboration [39| used Nf = 2 + 1 Clover fermions with a single mild stout smearing and a 
lattice spacing a = 0.076(2) fm. Finally, the LPHC collaboration @ obtained results using a hybrid action of domain 
wall valence quarks on a staggered sea on lattice of spatial length 2.5 fm and 3.5 fm at lattice spacing a = 0.124 fm. 

In Fig. [15] we compare our results on the strange octet baryons with those from the BMW, the PACS-CS, the 
QCDSF-UKQCD and the LHPC collaborations. Our results and the results by the PACS-CS and LHPC are not 
continuum extrapolated. The BMW results are extrapolated to the continuum limit and have larger errors than the 
rest. Nevertheless, there is an overall agreement, indicating that cut-off effects are small. In Fig. [TBI we compare our 
results on the strange decuplet baryons with the ones by PACS-CS and LHPC. Lattice results are in agreement except 
for the case of the S* where our results are consistently lower. Given the agreement of our results in the case of the 
51 this deviation cannot come from the mismatch in the strange quark mass. It is not clear what is the origin of this 
deviation for the S*, but the fact that the value obtained by PACS-CS at almost physical pion mass is high than the 
experimental value may indicate that the strange quark mass is larger than physical. In Fig. [T7] we show the masses 
for the strange baryons after extrapolating to the physical pion mass. Fort the results of this work we plot the values 
extracted using NLO HBxPT. Error shown on the twisted mass results is the estimate of the systematic error due to 
the chiral extrapolation, whereas the statistical errors are equal to the size of the symbols and are not shown. As can 
be seen, our results are in agreement with experiment except for f2, which is higher by 2%, just like the value found 
by PACS-CS. 

We also compare in Fig.[TS]our results for the charm baryons to those obtained using dynamical gauge configurations. 
All previous lattice computations of the mass of charm baryons used gauge confi gura tions produced with staggered 
sea quarks with a number of different actions for the valence quarks. In Refs. [4l|, [42| a Clover charm valence quark 
was used on MILC Nf = 2 + 1 gauge configurations at three values of the lattice spacing, a = 0.09,0.12,0.15 fm. In 
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FIG. 16: Comparison of the results of this work for the strange decuplet baryons to those obtained using domain wall valence 
quarks on a staggered sea 0].The notation is the same as that in Fig. [15] 
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FIG. 17: Lattice data in the strange quark sector extrapolated to the physical pion mass and the experimental values shown 
by the horizontal lines. For the twisted mass results of this work the chiral extrapolation was carried out using NLO HBxPT. 
We include in addition the results by the QCDSF-UKQCD collaboration [3i| (light blue crosses) . The rest of the notation is 
the same as in Fig. 1151 
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FIG. 18: Masses for charm baryons with spin 1/2 (left) and spin 3/2 (right) computed within lattice QCD and the experimental 
values shown by the horizontal lines. Our results for A c , S c , E cc , £*, H* c and Q ccc are extrapolated to the physical pion mass 
(shown with the red filled circles) whereas for the rest we give the results obtained at m-n = 260 MeV and /3 — 4.2 (shown 
with open red circles). We include results obtained usin g a number of hybrid actions with staggered sea quarks from Refs. 
(open blue squares), [4l], [Z^] (open green triangles) and [43T| (filled magenta triangles). 



Ref. [43[ Nf = 2 + 1 + 1 gauge configurations produced using the highly improved staggered quark (HISQ) action. 
The valence light-quark (up, down and strange) propagators arc generated using the clover impoved Wilson action. In 
order to reduce discretization artifacts a relativistic heavy-quark action was adopted for the charm quark. Finally in 
Ref. [40| domain wall fermions are used for the up, the down and the strange quarks on Nf = 2 + 1 improved Kogut- 
Susskind sea quarks at one value of the lattice spacing a = 0.12 fm. The relativistic Fermilab action was employed 
for the charm quark. We show the comparison of our results to those obtained in the aforementioned references in 
Fig. 1181 Our results for A c , S c , S cc , £*, 5* c and f2 ccc are extrapolated to the physical pion mass using a polynomial 
fit with up mi^-terms. For the spin 1/2 H c , S' c , f2 c and f2 cc we show the results obtained at m n = 260 MeV, the 
smallest value of the pion mass considered in this work on the lattice with the smallest lattice spacing at j3 = 4.2, for 
which cut-off effects arc smallest. As can be seen, our results are in agreement with the results of the other studies 
except for the S cc and with the experimental values. Although for S cc we find a value consistent with the result of 
the SELEX experiment, one has to study the pion mass dependence in order to reach a final conclusion. In Fig. [TH] 
we compare results for the spin 3/2 charm baryons. Our results for S*, f2* and f2* c are obtained at = 260 MeV 
at p = 4.2. There is good agreement among lattice results and with the known experimental values for £*, 2* and 
fi*. Thus the lattice results can be taken as a prediction for the masses of the other charm spin 3/2 baryons shown 
in the figure. 



VI. SUMMARY AND CONCLUSIONS 



In this work we have computed the strange and charm baryon masses using Nf = 2 twisted mass fermions. For the 
strange and charm sector we use an Osterwaldcr-Seiler valence quarks. The bare strange and charm valence quark 
mass is tuned by requiring that the physical values of the mass of the kaon and D-meson are reproduced after the 
lattice results are extrapolated at the physical value of the pion mass. 

We analyze gauge configurations for three values of the lattice spacings at the largest and smallest pion mass used 
in this study. Wc find that cut-off effects are small even in case of the charm baryons. This is a somewhat surprising 
result given that the Compton wave length of the D-meson mass is same order as the lattice spacing. 

Using simulations on two different volumes we obtained results that are consistent showing that any volume effects 
are smaller than our statistical accuracy. 

Another artifact of our lattice formulation is isospin breaking at finite lattice spacing. We have found that isospin 
breaking decreases with the lattice spacing and it is consistent with zero for a = 0.056 fm confirming the expected 
restoration of isospin symmetry. 

Our results on the strange quark sector are consistent with recent results using Clover improved fermions and 
domain wall fermions on a staggered sea. There is an overall agreement also in thee case of charm sector where we 
compare our results to other studies that used staggered sea quarks. The overall consistence among lattice results, 
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despite the different discretizations used, provides a strong validation of lattice QCD computations. Our results on 
the charm baryons reproduce the experimentally known values and thus provide an estimate for the mass of the n cc , 
Z* cc , and ^ccc- 
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VII. APPENDIX 



In Tablcs lXllXII and lXIII the strange octet baryon masses are collected for /3 = 3.9, j3 = 4.05 and /3 = 4.2, respectively. 
In Tables IXIII1 IXIVI and IXVI the strange dccuplet baryon masses for the three beta are collected. In Tables IXVI1 
IXVIII and IXVIIII we show the charm octet baryon masses for the three beta, while in Tables IXIXI IXXI and IXXII we 
show the charm decuplet baryon masses. In Table IXXIII we show the masses for extra spin-1/2 and spin-3/2 charm 
baryons at (3 = 4.2 at the lightest pion mass we have. 



L/a afih a[ii am s o 



am E + 



arriA 



32 
24 
.32 
24 
24 
24 
24 
24 
24 
32 
32 
24 
24 



0.0217 

0.0217 

0.0217 

0.0217 

0.0217 

0.0217 

0.015 

0.015 

0.015 

0.025 

0.030 

0.030 

0.030 



0.0030 
0.0040 
0.0040 
0.0064 
0.0085 
0.0100 
0.0064 
0.0085 
0.0100 
0.0040 
0.0040 
0.0064 
0.0085 



0.597(3 
0.610(5 
0.616(4 
0.628(5 
0.649(3 
0.654(4 
0.596(9 
0.641(5 
0.635(7 
0.623(5 
0.645(4 
0.654(5 
0.688(5 



0.577(4) 

0.600(3) 

0.596(5) 

0.610(6) 

0.631(2) 

0.643(5) 

0.585(6) 

0.598(11) 

0.627(8) 

0.606(5) 

0.630(5) 

0.640(6) 

0.651(9) 



0.607(7) 

0.622(4) 

0.628(7) 

0.640(9) 

0.668(4) 

0.666(6) 

0.633(9) 

0.663(6) 

0.645(15) 

0.633(9) 

0.664(6) 

0.674(7) 

0.703(6) 



0.657(4) 

0.671(3) 

0.675(4) 

0.687(5) 

0.697(3) 

0.696(5) 

0.658(11) 

0.674(5) 

0.659(14) 

0.687(6) 

0.729(6) 

0.732(5) 

0.739(6) 



0.628(3) 

0.638(2) 

0.635(3) 

0.650(4) 

0.659(2) 

0.668(3) 

0.597(12) 

0.621(6) 

0.622(7) 

0.651(4) 

0.688(4) 

0.682(7) 

0.706(4) 



0.561(3 
0.577(4 
0.574(4 
0.602(4 
0.619(3 
0.633(4 
0.588(4 
0.594(9 
0.614(6 
0.587(8 
0.606(5 
0.624(4 
0.626(9 



TABLE X: Strange octet baryon masses at /? = 3.9. 



L/a a^, h 



32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 



0.0178 

0.0178 

0.0178 

0.014 

0.014 

0.014 

0.0166 

0.0166 

0.0166 

0.019 

0.020 

0.020 

0.020 

0.025 



0.0030 
0.0060 
0.0080 
0.0030 
0.0060 
0.0080 
0.0030 
0.0060 
0.0080 
0.0060 
0.0030 
0.0060 
0.0080 
0.0060 



0.477(5 
0.500(5 
0.512(4 
0.463(5 
0.495(4 
0.494(8 
0.477(5 
0.496(5 
0.510(3 
0.501(6 
0.486(6 
0.510(6 
0.516(4 
0.522(5 



0.470(8 
0.483(5 
0.506(5 
0.456(9 
0.473(6 
0.508(9 
0.465(5 
0.483(5 
0.502(3 
0.495(8 
0.478(8 
0.493(5 
0.510(7 
0.501(6 



0.488(9) 

0.503(8) 

0.522(4) 

0.474(9) 

0.497(8) 

0.508(5) 

0.487(6) 

0.502(5) 

0.518(4) 

0.519(13) 

0.495(8) 

0.515(7) 

0.514(7) 

0.531(7) 



0.519(6 
0.529(5 
0.541(4 
0.496(6 
0.513(6 
0.512(7 
0.515(4 
0.526(4 
0.535(3 
0.532(6 
0.531(6 
0.544(6 
0.547(4 
0.572(5 



0.504(4 
0.519(4 
0.527(3 
0.482(4 
0.500(4 
0.507(3 
0.498(3 
0.511(4 
0.521(3 
0.524(6 
0.517(4 
0.534(4 
0.533(4 
0.556(4 



am,A 

0.455(5) 

0.482(5) 

0.496(5) 

0.445(6) 

0.479(4) 

0.485(6) 

0.451(4) 

0.480(5) 

0.496(3) 

0.483(5) 

0.461(5) 

0.492(5) 

0.500(5) 

0.501(5) 



TABLE XI: Strange octet baryon masses at /3 = 4.05. 



L/a a/j,h a/j,i am E o am E _ am E + am s o am s - am^ 

32 0.012 0.0065 0.402(4) 0.396(4) 0.405(5) 0.418(4) 0.406(4) 0.392(4) 

48 0.012 0.0020 0.362(3) 0.360(3) 0.368(4) 0.390(2) 0.381(2) 0.344(3) 

32 0.013 0.0065 0.406(7) 0.405(7) 0.417(7) 0.430(6) 0.416(6) 0.395(7) 

32 0.015 0.0065 0.411(4) 0.409(4) 0.413(4) 0.431(4) 0.421(4) 0.397(4) 

48 0.015 0.0020 0.374(2) 0.371(3) 0.380(3) 0.404(2) 0.397(2) 0.352(2) 

32 0.016 0.0065 0.417(7) 0.409(4) 0.413(4) 0.442(5) 0.421(4) 0.397(4) 



TABLE XII: Strange octet baryon masses at f3 — 4.2. 
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L/a afih afii 



am E »o 



amn 



32 
32 
24 
24 
24 
24 
24 
24 
24 
32 
32 
24 
24 



0.0217 

0.0217 

0.0217 

0.0217 

0.0217 

0.0217 

0.015 

0.015 

0.015 

0.025 

0.030 

0.030 

0.030 



0.0030 
0.0040 
0.0040 
0.0064 
0.0085 
0.0100 
0.0064 
0.0085 
0.0100 
0.0040 
0.0040 
0.0064 
0.0085 



0.700(7) 

0.726(4) 

0.716(12) 

0.749(7) 

0.759(5) 

0.773(8) 

0.737(24) 

0.751(10) 

0.730(18) 

0.730(9) 

0.724(11) 

0.772(6) 

0.782(11) 



0.713(5) 

0.718(6) 

0.726(14) 

0.749(7) 

0.755(4) 

0.772(8) 

0.736(7) 

0.736(20) 

0.757(19) 

0.723(12) 

0.743(10) 

0.769(6) 

0.736(7) 



0.700(6) 

0.715(6) 

0.717(17) 

0.740(11) 

0.755(5) 

0.773(9) 

0.746(7) 

0.755(7) 

0.746(16) 

0.710(14) 

0.739(9) 

0.760(11) 

0.746(7) 



0.710(6) 

0.739(4) 

0.737(8) 

0.754(8) 

0.767(3) 

0.784(5) 

0.728(15) 

0.723(11) 

0.720(15) 

0.747(12) 

0.782(7) 

0.797(7) 

0.812(7) 



0.702(5) 

0.721(4) 

0.714(8) 

0.734(8) 

0.738(5) 

0.749(7) 

0.708(10) 

0.659(18) 

0.667(10) 

0.738(8) 

0.765(8) 

0.781(6) 

0.751(11) 



0.798(3) 

0.809(3) 

0.805(7) 

0.815(6) 

0.820(3) 

0.834(4) 

0.772(16) 

0.771(9) 

0.771(11) 

0.827(7) 

0.870(6) 

0.874(6) 

0.771(9) 



TABLE XIII: Strange decuplet baryon masses at j3 = 3.9. 



L/a afj, h a/j,i 



amn 



32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 



0.0178 

0.0178 

0.0178 

0.014 

0.014 

0.014 

0.0166 

0.0166 

0.0166 

0.019 

0.020 

0.020 

0.020 

0.025 



0.0030 
0.0060 
0.0080 
0.0030 
0.0060 
0.0080 
0.0030 
0.0060 
0.0080 
0.0060 
0.0030 
0.0060 
0.0080 
0.0060 



0.597(11) 

0.593(8) 

0.606(6) 

0.589(12) 

0.586(11) 

0.588(10) 

0.570(12) 

0.582(12) 

0.615(4) 

0.592(16) 

0.611(7) 

0.599(10) 

0.610(6) 

0.598(11) 



0.589(10) 

0.586(8) 

0.607(7) 

0.595(8) 

0.574(13) 

0.606(7) 

0.578(7) 

0.578(10) 

0.609(6) 

0.576(13) 

0.605(8) 

0.591(9) 

0.619(6) 

0.589(12) 



0.590(13) 

0.582(8) 

0.606(6) 

0.603(6) 

0.582(12) 

0.592(6) 

0.561(8) 

0.568(14) 

0.604(7) 

0.551(19) 

0.610(7) 

0.595(11) 

0.605(6) 

0.577(15) 



0.596(8) 

0.582(7) 

0.609(6) 

0.578(8) 

0.568(10) 

0.574(9) 

0.573(6) 

0.579(6) 

0.601(5) 

0.579(8) 

0.606(7) 

0.594(11) 

0.613(8) 

0.613(7) 



0.588(6) 

0.583(7) 

0.598(5) 

0.570(6) 

0.559(11) 

0.574(6) 

0.567(6) 

0.577(7) 

0.586(5) 

0.583(9) 

0.596(8) 

0.591(9) 

0.597(6) 

0.613(8) 



0.661(7 
0.638(6 
0.651(6 
0.627(8 
0.623(7 
0.614(8 
0.630(5 
0.622(8 
0.648(4 
0.648(5 
0.671(5 
0.657(6 
0.660(5 
0.674(8 



TABLE XIV: Strange decuplet baryon masses at /? = 4.05. 



L/a ajXK apu 



amn 



32 0.012 0.0065 0.480(7) 0.482(6) 0.476(8) 0.468(6) 0.464(6) 0.498(6) 

48 0.012 0.0020 0.453(4) 0.451(4) 0.451(4) 0.444(3) 0.440(3) 0.487(3) 

32 0.013 0.0065 0.501(11) 0.500(10) 0.495(14) 0.494(9) 0.489(9) 0.522(10) 

32 0.015 0.0065 0.478(6) 0.487(5) 0.477(6) 0.478(5) 0.475(5) 0.513(4) 

48 0.015 0.0020 0.458(4) 0.455(3) 0.452(2) 0.440(4) 0.448(4) 0.505(3) 

32 0.016 0.0065 0.507(11) 0.487(5) 0.503(12) 0.504(10) 0.505(9) 0.540(10) 



TABLE XV: Strange decuplet baryon masses at f3 = 4.2. 



L/a 


a\ih 


a[ii 


am E + 


om £ o 


am E ++ 


am„++ 




am A + 


24 


0.240 


0.0040 


1.100(6) 


1.105(13) 


1.102(14) 


1.532(5) 


1.528(4) 


1.015(8) 


24 


0.240 


0.0064 


1.117(6) 


1.059(20) 


1.122(6) 


1.552(3) 


1.533(4) 


1.045(5) 


24 


0.240 


0.0085 


1.131(6) 


1.125(5) 


1.135(5) 


1.555(4) 


1.541(5) 


1.063(3) 


24 


0.240 


0.0100 


1.139(4) 


1.131(5) 


1.138(5) 


1.559(3) 


1.551(2) 


1.070(3) 


24 


0.250 


0.0040 


1.128(7) 


1.120(9) 


1.099(14) 


1.591(5) 


1.575(5) 


1.055(6) 


24 


0.270 


0.0040 


1.154(5) 


1.157(5) 


1.161(6) 


1.640(4) 


1.628(3) 


1.069(4) 


24 


0.270 


0.0064 


1.172(7) 


1.164(5) 


1.176(7) 


1.648(5) 


1.634(4) 


1.096(5) 


24 


0.270 


0.0085 


1.192(4) 


1.181(4) 


1.188(5) 


1.658(3) 


1.651(3) 


1.115(3) 


24 


0.270 


0.0100 


1.196(4) 


1.189(4) 


1.197(4) 


1.660(3) 


1.649(25) 


1.122(3) 


24 


0.300 


0.0040 


1.203(9) 


1.214(5) 


1.221(5) 


1.736(6) 


1.731(4) 


1.135(4) 


24 


0.300 


0.0064 


1.227(4) 


1.162(22) 


1.236(5) 


1.747(4) 


1.722(6) 


1.148(5) 


24 


0.300 


0.0085 


1.244(4) 


1.236(4) 


1.244(4) 


1.759(3) 


1.736(6) 


1.169(3) 


24 


0.300 


0.0100 


1.251(3) 


1.237(5) 


1.249(4) 


1.762(2) 


1.751(3) 


1.172(3) 



TABLE XVI: Charm spin-1/2 baryon masses at p = 3.9. 
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L/a au,h dUi am^+ am s o am„++ am^++ am~+ 



32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 



0.170 
0.170 
0.170 
0.200 
0.200 
0.200 
0.210 
0.210 
0.210 
0.230 
0.230 
0.230 
0.260 
0.260 
0.260 



0.0030 
0.0060 
0.0080 
0.0030 
0.0060 
0.0080 
0.0030 
0.0060 
0.0080 
0.0030 
0.0060 
0.0080 
0.0030 
0.0060 
0.0080 



0.843(4 
0.844(4 
0.852(4 
0.900(4 
0.903(4 
0.905(5 
0.917(4 
0.919(4 
0.926(3 
0.952(5 
0.952(4 
0.959(5 
1.004(5 
1.001(5 
1.010(6 



0.841(5) 

0.831(9) 

0.852(5) 

0.888(7) 

0.886(10) 

0.907(5) 

0.912(6) 

0.911(4) 

0.925(4) 

0.942(7) 

0.940(11) 

0.960(5) 

1.005(10) 

0.992(11) 

1.011(7) 



0.839(6 
0.836(8 
0.852(7 
0.894(6 
0.891(8 
0.908(8 
0.908(7 
0.921(4 
0.926(4 
0.948(6 
0.945(8 
0.959(5 
1.003(5 
0.997(8 
1.012(8 



1.134(5 
1.149(3 
1.146(2 
1.244(3 
1.255(2 
1.250(2 
1.282(3 
1.286(3 
1.288(2 
1.345(3 
1.355(3 
1.349(3 
1.444(3 
1.457(3 
1.449(3 



1.137(4 
1.139(4 
1.142(3 
1.244(5 
1.250(3 
1.247(3 
1.276(4 
1.279(3 
1.283(2 
1.342(5 
1.353(3 
1.349(3 
1.441(5 
1.451(4 
1.446(5 



0.774(5 
0.791(5 
0.791(6 
0.828(6 
0.843(5 
0.848(5 
0.845(6 
0.859(5 
0.865(4 
0.883(5 
0.895(5 
0.899(5 
0.935(5 
0.946(5 
0.955(3 



TABLE XVII: Charm spin-1/2 baryon masses at /3 = 4.05. 



L/a afih ai±i am s + am E o am E ++ am„++ am„+ am A + 

32 0.130 0.0065 0.696(6) 0.694^6) 0.698(6) 0.932(4) 0.927(5) 0.653(4) 

32 0.160 0.0065 0.733(6) 0.731(6) 0.734(7) 0.999(4) 0.997(5) 0.688(4) 

32 0.185 0.0065 0.778(7) 0.776(6) 0.779(7) 1.085(4) 1.082(5) 0.731(4) 

32 0.210 0.0065 0.821(7) 0.819(7) 0.822(7) 1.168(4) 1.156(5) 0.774(4) 

48 0.136 0.0020 0.653(3) 0.656(3) 0.652(3) 0.899(2) 0.898(2) 0.603(2) 

48 0.170 0.0020 0.716(3) 0.719(3) 0.715(4) 1.017(2) 1.016(2) 0.663(3) 



TABLE XVIII: Charm spin-1/2 baryon masses at P = 4.2. 



L/a a[ih afj,i am E »+ am^*o 



++ am-z 



amn a 



24 
24 
24 
24 
24 
24 
24 
24 
24 
24 
24 
24 
24 



0.240 
0.240 
0.240 
0.240 
0.250 
0.270 
0.270 
0.270 
0.270 
0.300 
0.300 
0.300 
0.300 



0.0040 
0.0064 
0.0085 
0.0100 
0.0040 
0.0040 
0.0064 
0.0085 
0.0100 
0.0040 
0.0064 
0.0085 
0.0100 



1.148(10) 

1.159(10) 

1.175(8) 

1.184(6) 

1.173(8) 

1.201(8) 

1.209(10) 

1.231(4) 

1.239(4) 

1.263(5) 

1.267(5) 

1.281(4) 

1.292(3) 



1.142(16) 

1.151(9) 

1.164(8) 

1.184(4) 

1.182(9) 

1.204(7) 

1.210(5) 

1.224(5) 

1.234(4) 

1.245(15) 

1.260(5) 

1.274(5) 

1.283(4) 



1.147(15) 

1.166(11) 

1.174(10) 

1.181(6) 

1.173(11) 

1.210(5) 

1.225(6) 

1.230(4) 

1.238(5) 

1.260(6) 

1.275(6) 

1.280(4) 

1.282(6) 



1.572(6 
1.580(6 
1.594(5 
1.599(3 
1.602(7 
1.671(7 
1.680(5 
1.692(5 
1.694(4 
1.775(6 
1.779(5 
1.794(3 
1.795(3 



1.564(6 
1.572(5 
1.578(7 
1.591(3 
1.606(8 
1.668(5 
1.672(5 
1.686(3 
1.690(4 
1.764(5 
1.768(5 
1.770(8 
1.789(3 



1.989(3) 
1.991(4) 
1.997(4) 
1.999(3) 
2.043(4) 
2.130(3) 
2.133(3) 
2.136(4) 
2.141(3) 
2.269(3) 
2.270(3) 
2.274(4) 
2.277(2) 



TABLE XIX: Charm spin-3/2 baryon masses at /3 = 3.9. 
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L/a afih afii 



am- S t+ onij.o 



«++ am=* ++ om = . + amo 



32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 
32 



0.170 
0.170 
0.170 
0.200 
0.200 
0.200 
0.210 
0.210 
0.210 
0.230 
0.230 
0.230 
0.260 
0.260 
0.260 



0.0030 
0.0060 
0.0080 
0.0030 
0.0060 
0.0080 
0.0030 
0.0060 
0.0080 
0.0030 
0.0060 
0.0080 
0.0030 
0.0060 
0.0080 



0.889(6 
0.887(4 
0.885(5 
0.940(6 
0.928(6 
0.936(5 
0.956(6 
0.955(5 
0.957(4 
0.989(6 
0.984(4 
0.986(5 
1.039(6 
1.033(5 
1.034(5 



0.889(6 
0.883(5 
0.890(7 
0.940(6 
0.937(4 
0.945(6 
0.956(6 
0.954(4 
0.962(4 
0.989(6 
0.987(4 
0.995(6 
1.038(6 
1.032(5 
1.039(6 



0.890(4) 

0.868(12) 

0.884(5) 

0.941(4) 

0.918(11) 

0.935(5) 

0.958(4) 

0.948(5) 

0.958(4) 

0.992(4) 

0.967(11) 

0.986(5) 

1.036(7) 

1.016(11) 

1.032(10) 



1.170(6 
1.184(3 
1.176(3 
1.277(4 
1.288(3 
1.279(3 
1.310(4 
1.316(3 
1.316(3 
1.376(4 
1.387(3 
1.378(3 
1.473(4 
1.485(3 
1.472(3 



) 1.174(5) 


1.468(3) 


) 1.172(6) 


1.476(2) 


) 1.175(4) 


1.468(3) 


) 1.281(4) 


1.616(3) 


) 1.281(3) 


1.623(2) 


) 1.281(3) 


1.612(3) 


) 1.314(4) 


1.663(3) 


) 1.310(4) 


1.667(2) 


) 1.313(3) 


1.663(3) 


) 1.379(4) 


1.761(3) 


) 1.380(4) 


1.768(2) 


) 1.380(3) 


1.758(3) 


) 1.476(4) 


1.905(3) 


) 1.478(4) 


1.910(3) 


) 1.472(5) 


1.898(4) 



TABLE XX: Charm spin-3/2 baryon masses at j3 = 4.05. 



L/a aiih aui am T *+ am^.o am T *++ arriw* ++ am=. + amn rr . r 



32 0.130 0.0065 0.730(8) 0.730(8) 0.727(8) 0.958(6) 0.963(6) 1.191(3) 

32 0.160 0.0065 0.763(8) 0.763(8) 0.761(8) 1.025(6) 1.029(6) 1.287(3) 

32 0.185 0.0065 0.805(8) 0.804(8) 0.802(8) 1.107(6) 1.109(6) 1.408(3) 

32 0.210 0.0065 0.845(9) 0.845(8) 0.843(8) 1.188(6) 1.193(6) 1.526(3) 

48 0.136 0.0020 0.686(4) 0.688(3) 0.683(4) 0.925(2) 0.926(2) 1.166(1) 

48 0.170 0.0020 0.744(4) 0.746(3) 0.741(4) 1.039(2) 1.039(2) 1.333(1) 



TABLE XXI: Charm spin-3/2 baryon masses at p = 4.2. 



L/a au s au c aui am-= a am-^i amn c amn al , ara-s.* amn* amn* 
48 0.015 0.17 0.0020 0.708(2) 0.745(3) 0.771(2) 1.044(1) 0.770(3) 0.794(2) 1.065 C (2) 



TABLE XXII: Strange-charm spin-1/2 and 3/2 baryon masses at /3 — 4.2 at the tuned heavy quark masses. 



